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ABSTRACT 

A  theory  for  predicting  the  stress-strain  characteristics 
of  polymeric  solid3  is  developed  in  terms  of  a  description  of 
microdefect  formation.  The  process  of  irreversible  change  in 
these  solids  is  assumed  to  be  a  combination  of  nucleation  of  sub- 
microscopic  defects  at  stress  inhomogeneities  and  their  subsequent 
growth  to  macroscopic  dimensions.  Straining  results  in  the 
simultaneous  generation  of  free  volume  and  macrodefects  so  that 
catastrophic  failure  can  occur  through  either  a  general  yielding 
of  the  material  or  by  brittle  fracture. 

It  is  assumed  that  nucleation  of  sub-microscopic  defects  is 
an  activated  process  and  that  defect  growth  is  one  dimensional 
and  linear.  The  total  strain  is  expressed  as  the  sum  of  an 
elastic  recoverable  strain  and  a  nonlinear,  nonrecoverable  strain 
and  expressions  are  obtained  for  the  stress  as  a  function  of  time, 


*On  leave  from  the  Laboratorio  di  Ricerca  su  Tecnologia  dei 
Polimeri  e  Reologia  del  C.N.R.  Napoli,  Italy. 


temperature  and  loading  history.  The  criterion  for  yielding  is 
defined  in  terms  of  free  volume  concepts  and  a  criterion  for 
brittle  failure  is  defined  in  terms  of  a  critical  defect  size. 

The  parameters  of  the  resulting  model  are  calculated  for  poly¬ 
phenylene  oxide  polymers  based  on  constant  rate  of  loading  experi¬ 
ments  and  then  the  general  creep  behavior,  including  the  time 
required  under  constant  load  for  cold  flow,  is  predicted.  Experi¬ 
mental  data  are  shown  to  agree  with  these  predictions. 
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Introduction 

At  temperatures  well  below  the  primary  glass  transition, 
most  organic  polymers  exhibit  either  brittle  or  ductile  failure, 
depending  upon  the  load  and  temperature  history  imposed  on  the 
material.  When  ductile,  these  polymers  are  tough  and  resistant 
to  impact  and  when  brittle  they  are  not.  Polyblending  with  finely 
dispersed  rubber  particles  often  results  in  large  increases  in 
toughness,  while  polyblending  with  coarse  rigid  particles  often 

results  in  increased  brittleness . 

The  difficulty  in  describing  the  stress-strain  behavior 
complicated  by  many  factors,  of  which  the  formation  and  growth 
of  defects  during  loading  is  perhaps  the  most  important.  It  is 
well  known  that  the  presence  of  stress  inhomogeneities  and/or 
finely  dispersed  second  phases  can  induce  crazes,  cold  flow  and 
a  multitude  of  interacting  cracks.  It  is  the  purpose  of  this  paper 
to  offer  a  rational  theory  to  describe  the  stress-strain  behavior 
of  glassy  polymeric  solids. 

Theory 

The  process  of  irreversible  change  in  a  polymeric  solid  is 
assumed  to  be  a  combinatJ ''n  of  nucleation  of  sub-microscopic 
defects  at  stress  inhomogeneities  and  their  subsequent  growth  to 

These  defec  might  be  either  true  cracks, 


macroscopic  dimensions. 


in  which  case  continued  straining  of  the  material  results  in  the 
generation  of  a  volume  change  in  the  solid  (1)  or  they  might  be 
crazes,  in  which  case  continued  straining  of  the  material  results 
in  the  generation  of  free  volume  (2-4).  In  any  case,  they  lead 
to  catastrophic  failure  through  either  a  general  yielding  of  the 
material  or  a  brittle  fracture. 

In  developing  any  theory  of  failure  it  is  necessary  to  set 
criteria  for  both  types  of  failure  e*nd  then  predict  which  will 
occur  first  under  a  given  set  of  loading  conditions. 

In  this  paper  we  deal  specifically  with  amorphous  polymeric 
solids  which  tend  to  craze  under  the  influence  of  an  external 
load.  It  has  been  observed  that  under  a  unidirectional  tensile 
load  crazes  will  nucleate  and  grow  perpendicular  to  the  direction 
of  loading.  Under  certain  conditions  they  tend  to  overlap  and 
coalesce  to  form  a  nearly  continuous  zone  of  crazed  material 
which  ultimately  leads  to  a  general  cold  drawing,  or  macroscopic 
yielding,  of  the  solid.  This  condition  can  be  observed  experi¬ 
mentally  by  the  development  of  necking  and  optical  birefringence 
in  the  solid,  or  less  precisely,  by  the  appearance  of  an  apparent 
maximum  in  the  engineering  stress-strain  curve.  Under  certain 
conditions, these  crazes  can  grow  to  a  greater  length  than  can  be 
tolerated  by  the  material  and  the  craze  rapidly  changes  to  a 
macroscopic  crack  which  propagates  catastrophically,  causing 
failure  in  a  brittle  manner. 

It  is  assumed  that  the  nucleation  of  craze  sites  is  an 
activated  process.  Based  on  the  work  of  Zhurkov  et  al.  (5), 


one  can  empirically  justify  the  following  form  for  the  initial 
nucleation  rate: 

dN  /  AEn\  i  i  ,  )  (D 

§£  “  No  ^exp  -  -55  J  sinh\  c*  / 

where  is  a  constant,  AE„  is  an  activation  energy  and  o*/(n-l) 
is  a  characteristic  stress  constant.  For  very  long  ter*  loading 
•Xhurkov's  data  show  that  it  is  necessary  to  account  for  an 
exponential  decrease  in  the  nucleation  rate  using  a  decay  factor 
exp  -  t/X  (o,  T) .  However,  Equation  (1)  is  satisfactory  for  short 

term  stress-strain  behavior . 

Once  craze  sites  have  been  nucleated,  their  growth  can  be 
represented  by  a  one  dimensional,  linear  growth  process  (6).  It 
is  assumed  that  the  growth  rate  is  given  by: 
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where  G  is  a  constant,  AE-  an  activation  energy  and  o*  a 
o  ° 

characteristic  stress  constant. 

The  isothermal  rate  of  formation  of  crate  volume  (dVf/dt)T 
can  then  be  obtained  from  the  product  of  nucleation  and  growth 


rates  as: 
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where  B(T)  =  $q  Nq  Gq  exp  -  - ^ - —  (4) 

and  $Q  is  a  proportionality  constant. 

In  the  glassy  state,  one  can  neglect  pure  viscous  flow  and 
following  the  suggestion  of  DiBenedetto  and  Trachte  (7)  ,  the 
total  strain  cm  be  expressed  as  the  sum  of  two  terms,  a 
linear  elastic  (and  recoverable)  strain  and  a  non-linear 
(and  irrecoverable)  strain.  The  elastic  strain  is  equated 
to  the  initial  elasticity  of  the  material  a/ Eg  and  the  irrecover 
able  strain  is  equated  to  the  strain  caused  by  the  additional 
craze  volume: 


£  -  es  +  ec 

0  V£(t) 

■  EJ  +  TT^w T  (5) 

where  a  is  the  stress,  Eg  the  initial  elastic  modulus,  u  the 
Poisson  ratio  and  (t)  the  craze  volume  obtained  from 
integration  of  equation  (3)  . 

Equation  (5)  gives  the  relationship  between  stress, 
strain,  temperature  and  time,  in  order  to  fix  a  point  of  failure 
for  the  material,  failure  criteria  must  be  established.  It 
is  assumed,  in  accordance  with  free  volume  concepts  (2,3), 
that  the  criterion  for  yielding,  or  the  onset  of  cold  flow, 


is  given  by: 


Vf(ty)  =  (Aa)  <Tm-T)  (6) 

where  (Aa)  is  a  constant  related  to  the  difference  in  thermal 
expansion  coefficients  of  the  glassy  and  fluid  states  and  Tm 
is  the  maximum  temperature  at  which  craze  sites  can  nucleate. 

The  temperature  Tm  has  been  related  to  the  glass  tempera¬ 
ture  Tg  (7)  and  sometimes  is  considered  to  be  slightly  higher 
than  T  (8) .  In  this  paper  it  is  left  as  a  measurable  parameter 
It  is  also  assumed  that  the  criterion  for  brittle 
failure  is  determined  by  a  critical  flaw  length  i£(T,  p^,) 
that  is  primarily  a  function  of  temperature  and  flaw 
density. 

Thus,  whether  the  material  yields  before  brittle  fracture 
depends  upon  whether  the  formation  of  crazes  develops  sufficient 
free  volume  Vf  (*-j)  before  reaching  a  critical  flaw  size 
Which  criterion  is  reached  first,  depends  upon  the  relative 
rates  of  nucleation  and  growth.  The  stress-strain  behavior 

under  some  well-defined  loading  conditions  will  be  considered 

!  - 

next. 

A.  Creep 

In  a  creep  experiment  the  stress  is  constant  at  c0.  The 
integration  of  equation  (3)  and  substitution  of  Vf (t)  into 
equation  (5)  gives: 
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The  first  term  on  the  right  side  of  equation  (7)  represents 
the  recoverable  elastic  deformation  and  the  second  term  the 
irrecoverable  deformation  associated  with  craze  formation.  The 
use  of  equation  (6)  for  the  yield  criterion  leads  to  the  follow¬ 
ing  expression  for  the  time  to  yield  ty  (i.e.,  the  time  required 


to  initiate  cold  drawing) 
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Integration  of  equation  (2)  gives  the  defect  size  at  the 
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yield  point  G(ty)  as: 

/  AE_\  /  ^o\ 

G(ty)  -  G0  (exp  --ef)  (sinh 

If  g (ty)  >  fcj  brittle  failure  occurs  before  cold  drawing. 
B.  Constant  Loading  Rats 

In  this  case  the  stress  is  a  linear  function  of  time, 
a  -  Es  rn  t,  and  the  rate  of  loading  or  is  constant  and  equal 
The  integration  of  equation  (3)  gives: 


tO  Eg  V 


B_(T) 

(n 


T)q*a  c_ 
i-l)a2 


no  (2-n)a 

cosh  cosh  0y— 


2n 


2  (2-n) 


-  cosh 


a 

"a* 


n(2-nT 


+  1 


From  experimental  observation  it  is  found  that  n  >  1  and 
a  /a*  »>  1  so  that  for  all  practical  purposes  over  most  of 
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the  stress-strain  curve, equation  10  reduces  to: 
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Substitution  of  equation  (11)  into  equation  (5)  gives: 


et(a) 


B (T)  a*2  (exp  2°  ) 
4  a2  n(n-l)  (l-2y) 


(12) 


Since  the  criterion  for  yielding  is  Vf(a y)  »  (Aa) (T^-T) ,  the 
yield  stress  can  be  obtained  from  equation  (11)  as: 


The  sensitivity  of  yield  stress  to  loading  rate  is  given 

by: 


d(Ty -  -  isl  d4) 

d  in  or  n 

which  predicts  that  yield  stress  is  a  linear  function  of  rate 
of  loading.  The  temperature  dependence  of  yield  stress  is 
given  by: 
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which  shows  that  the  temperature  dependence  of  yield  stress 
is  controlled  primarily  by  the  difference  in  activation 
energies  between  the  nucleation  and  growth  processes. 

The  defect  size  at  the  yield  point  can  be  obtained  from 


integration  of  equation  (2) : 
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C.  Constant  Strain  Rate 

The  maior  difficulty  in  analyzing  constant  strain  rate 
data  is  that  the  rate  of  formation  of  craze  volume  is  defined 
as  an  integral  of  stress  over  time  (equation  3)  and  since  the 
stress-strain  curve  is  non-linear,  the  dependence  of  stress 
on  time  is  not  known.  Cne  can,  however,  define  the  stress- 
strain  relation  as  a  ■  Eg(<j)(r0t)  where  Eg  (a)  is  a  stress 
dependent  modulus  and  rQ  is  the  rate  of  straining.  Prom 
observation  of  experimental  stress-strain  curves,  a  reasonable 
expression  for  the  time  dependence  of  the  stress  to  nearly 
the  yield  point  is: 


da  •  _  5a 

3t  ESo  ro  exP  -  (17) 

where  6  is  a  constant  of  the  order  of  0.01  to  0.1.  The  use  of 
equation  (17)  in  equation  (3)  leads  to 

.  B(t)  a*2  exp  <* 

V  s  a* 

f  4  Es<j'r  roMn+5-l)  (n-25)  (18) 

and 

day  .  2a* 

d  in  rQ  “  H+2T  (19) 

As  long  as  the  stress-strain  curve  is  nearly  linear, 
constant  rate  of  loading  and  constant  rate  of  straining  result 
in  nearly  the  same  behavior  with  the  yield  stress  being 
slightly  less  dependent  on  strain  rate  than  on  stress  rate. 


-s 
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Experimental  Results 

The  material  studied  was  General  Electric's  Grade  631-111 

polyphenylene  oxide  with  a  glass  transition  temperature  of 

T  *  210°C,  a  Poisson's  ratio  of  0.35  and  a  structural  formulas 
9 


Test  specimens  were  compression  molded  at  285°C  followed  by 

a  slow  cooling  to  room  temperature.  All  samples  were  annealed 

below  T  to  minimize  molding  stresses.  Standard  ASTM  tensile 
9 

tests  were  carried  out  on  a  modified  Instron  Testing  machine 
which  could  be  controlled  for  both  constant  load  and  constant 
rate  of  loading  as  well  as  constant  rate  of  deformation.  Tests 
were  conducted  at  temperatures  ranging  from  25°C  to  200°C  and 
at  strain  rates  ranging  from  0.00526  to  0.526  in/in-min. 

A  typical  stress-strain  curve  is  shown  in  Figure  1.  The 
yield  point  was  assumed  to  be  the  maximum  in  the  curve.  If 
failure  occurred  before  this  maximum,  it  was  defined  as 
brittle  failure.  The  clastic  strain  was  measured  from  the 
initial  elastic  modulus  and  the  irrecoverable  deformation  was 
associated  with  the  non-linear  portion  of  the  stress-strain 
curve.  A  summary  of  the  yield  stress  data  at  constant  loading 
rate  and  constant  deformation  rate  over  a  range  of  temperatures 
is  given  in  Figure  2.  Time-temperature  superpositioning  at  a 
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reference  temperature  of  34 °C  (9)  was  used  to  condense  the 
data  and  the  required  shift  factor  a^  is  given  by: 

log10  aT  =  ‘  28*° 

where  T  in  the  temperature  in  degrees  Kelvin. 

The  effect  of  temperature  on  yield  stress  is  shown  in 
Figure  3.  According  to  the  yield  criterion  (equation  6) , 
the  irreversible  strain  at  the  yeild  point  should  be  a  linear 
function  of  temperature  and  independent  of  loading  rate. 

Table  1  shows  that  it  was  independent  of  loading  rate  and 
its  temperature  dependence  can  be  expressed  as: 

(ep)y  *  V(V  (1’2y)  *  3*09  x  10~5  (525“T) 

in  agreement  with  theory.  Figure  4  illustrates  the  appearance 
of  the  crazes  for  a  variety  of  test  conditions.  An  average 
defect  size  was  measured  from  such  pictures  and  the  results 
are  shown  in  Figure  5.  Defect  size  is  seen  to  be  a  complex 
function  of  strain  rate  and  temperature.  As  long  as  the 
defect  size  was  below  the  locus  of  points  representing 
the  failure  was  by  yielding,  while  if  the  defect  size  was 
greater  than  l*  the  failure  was  brittle.  It  is  reasonable 
to  expect  that  the  critical  crack  size  is  a  function  of 
temperature,  defect  density  and  defect  size  distribution  so 
that  it  is  probably  more  realistic  to  imagine  a  narrow  zone 
about  the  locus  of  points  representing  l*  which  divides 
the  brittle  and  ductile  zones. 
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From  the  data  shown  in  Figures  2  to  5  and  using 
equations  4  and  12  to  16,  one  may  evaluate  all  of  the 
constants  required  to  define  the  stress-strain  behavior 
of  the  material.  These  values  are  listed  below: 

Aa  *  3.09  x  10"5  °K_1 

T  *  525°K 

m 

a*  **  334  psi 

n  »  2.3 

AEn  ■  41,700  cal/g-mole 

AEg  »  32,300  cal/g-mole 

8oGoNo  *  2  x  lo2°  -  100% 

The  values  for  the  activation  energies  are  reasonable 
for  the  assumed  processes  and  the  value  of  a*  is  very  low 
compared  to  the  yield  stresses  in  the  material.  Because  of  this 
latter  fact  the  constant  @0GQN0  is  susceptible  to  large  error. 
The  limits  of  +  100%  in  30GQN0  correspond  to  a  +  100  psi  error 
in  load  measurement  or  a  +  1°C  in  temperature  control. 

The  primary  value  of  such  a  model  is  to  predict  the 
mechanical  behavior  in  other  modes  of  loading,  based  on  the 
data  from  constant  loading  rate  experiments.  By  using 
equations  (7)  to  (9),  one  can  calculate  the  irreversible 
deformation,  the  time  to  yield  and  the  defect  size  at  the 
yield  point  under  a  constant  load.  Creep- recovery  experiments 
were  run  in  order  to  experimentally  determine  these  quanti¬ 
ties.  The  theoretical  curves  and  the  experimental  data  are 


14 


shown  in  Figures  6  to  8.  The  irrecoverable  deformation  as  a 
function  of  time  was  calculated  from  the  second  term  on  the 
right  hand  side  of  equation  (7).  The  experimental  values 
shown  on  Figu  e  6  agreed  with  theory  to  within  experimental 
accuracy.  Similarly,  the  experimental  values  for  time  to 
yield  agreed  with  those  calculated  from  equation  (8) ,  as 
shown  in  Figure  7.  When  the  failure  was  in  the  brittle  mode, 
it  occurred  earlier  than  the  calculated  value  for  yield  time. 
These  points  are  marked  by  the  symbol  B  in  Figure  7.  Equation  9 
is  used  to  calculate  average  defect  size  at  yield  and  the 
theoretical  curves  are  plotted  in  Figure  8.  The  values  of 
lC  shown  in  Figure  8  were  obtained  from  Figure  5.  As  predicted, 
when  the  defect  size  was  above  A*  a  brittle  failure  was 
encountered.  When  ductile  failure  was  attained,  the  calculated 
values  of  defect  size  agreed  with  experiment. 

Concluding  Remarks 

A  theory  for  predicting  the  stress-strain  characteristics 
of  polymeric  solids  has  been  presented.  The  process  of  irrevers¬ 
ible  change  is  assumed  to  be  a  combination  of  nucleation  of 
submicroscopic  defects  at  stress  inhomogeneities  and  their 
subsequent  growth  to  macroscopic  dimensions.  The  total  strain 
is  expressed  as  the  sum  of  an  elastic  recoverable  strain 
and  a  non-linear  non-recoverable  strain  and  expressions  are 
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obtained  for  the  stress  as  a  function  of  time,  temperature 
and  loading  history.  The  criterion  for  yielding  is  defined 
in  terms  of  free  volume  concepts  and  a  criterion  for  brittle 
failure  is  ^fined  in  terms  of  a  critical  defect  size. 

It  ha^  been  shown  that  the  resulting  model  quantitatively 
describes  the  mechanical  behavior  of  polyphenylene  oxide  and 
that  by  using  stress-strain  data  obtained  at  constant  rate  of 
loading,  one  can  predict  the  creep  bcnavior  of  the  material 
at  constant  load.  A  description  of  the  ductile-brittle 
transition  has  also  been  presented  and  it  has  been  shown  that 
one  can  qualitatively  define  the  locus  of  points  separating 
the  ductile  from  the  brittle  behavior. 

Concerning  this  latter  point,  it  is  well  known  that 
brittle  failure  is  a  stochastic  process  and  that  one  cannot 
predict  with  absolute  certainty  the  time  for  brittle  failure. 
This  implies  that  the  brittle-ductile  transition  cannot  be 
represented  by  a  single  locus  of  points,  but  rather  the  two 
regions  are  separated  by  a  narrow  transition  zone  within 
which  both  modes  of  failure  are  more  or  less  equally  probable. 
A  second  point  is  that  the  concept  of  a  simple  critical  crack 
dimension  for  brittle  failure  is  oversimplified  and  the 
criterion  for  brittle  failure  must  certainly  be  examined  in 
terms  of  stress  intensity  factors.  Both  of  these  points  are 
presently  under  study. 
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Table  1.  The  Irrecoverable  Strain  at  Yield,  ep  ■  Vf (Ty)/ (l-2y) 
as  a  Function  of  Temperature  and  Strain  Rate. 
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1.68 

1.68 

1.60 

115 

1.35 

1.40 

1.35 

1.52 

130 

1.30 

1.10 

1.18 

1.20 

1.22 

1.35 

1.25 

140 

1.18 

1.15 

1.13 

_ 1 

1.22 

e  in/in-min 


1 


■  -  •■•«««# 


LIST  OF  FIGURES 


r 

r 

r 

i 

\ 


i' 

|. 


Figure  1.  A  typical  engineering  stress-strain  curve  for  PPO. 
Figure  2.  Effect  of  strain  rate  on  yield  stress. 

Figure  3.  Effect  of  temperature  on  yield  stress. 

Figure  4.  Illustration  of  the  appearances  of  crazes. 

Figure  5.  Effect  of  tenqperature  and  strain  rate  on  average 


defect  size. 

Figure  (,  Non-recoverable  deformation  in  creep  as  a  function 
of  time. 

Figure  7.  Time  to  yield  under  constant  load. 

Figure  8.  The  effect  of  creep  stress  on  the  average  defect 
size. 
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